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ABSTRACT 

We present a systematic analysis of ground state properties of a tunneling 
system coupled to an Ohmic bath. The nonadiabatic effects induced by 
coupling with tunneling particle , on the ground state properties of phonons 
are taken into account. For this we have considered both displacement and 
deformation of phonons via a variational treatment. The symmetry breaking 
transition , or the delocalization-localization transition of tunneling particle 
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is presented and is compared with earlier known results. We also show that 
in the critical regime a gap will be opened near the zero momentum in the 
spectrum of the phonons. The modified phonon spectrum along with the 
phonon displacement field correlation functions have been calculated. It is 
shown that in the critical regime the tunneling system stabilizes the lattice 
against long wavelength fluctuations. 
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1. Introduction 



The statics and dynamics of two state tunneling system coupled to an en- 
vironment is of fundamental interest both in physics and chemistry [1]. This 
rather simple system has been used to study the motion of impurities (atoms 
or groups of atoms) in solids and glasses which play an important role in de- 
termining low temperature properties [2,3] , defect motion in alkali halides 
[4] , charged particle motion in metals, magnetic impurities in metals [5,6] 
and molecular transition in liquids [7] etc. In all examples mentioned above 
tunneling is a many-body problem in the sense that the considered parti- 
cle couples to a large number of degrees of freedom arising from medium or 
environment. In recent years there has been a growing intrest in the influ- 
ence of dissipation on the behaviour of macroscopic quantum variables and 
quantum tunneling out of metastable state. In particular, the trapped mag- 
netic flux threading a superconducting ring interrupted by weakly coupled 
Josephson junction is expected to show macroscopic quantum coherence. In 
this case the problem of tunneling of quantized magnetic flux reduces to 
quantum motion of a particle in a double well potential [1]. The existence 
of chiral molecules by a superselection rule which origanates from the ever 
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present coupling of the molecule to the radiation field can be understood 
in the frame work of two state system [8]. Here the two states are noth- 
ing but the left and right handed configurations of molecule. Recently the 
two state system formulation has been applied to an understanding of c-axis 
resistivity in high-T c layered oxides and transverse resistivity in quasi one- 
dimensional compounds [9,10] . Electrical transport along c-axis in layered 
oxides is pictured as a coherent interplanar tunneling between neighbouring 
layers blocked by repeated intraplanar incoherent (inelastic) scattering. Here 
two states correspond to two nearest neighbour layers. The ground state 
properties of two level system are also related to anisotropic Kondo problem 
[1] and one dimensional inverse square Ising model [11,12] . In particular 
the delocalization-localization transition in two state system is related to the 
existence of ferromagnetic phase transition in the one-dimensional inverse 
square Ising model and several rigorous results for this problem are known. 

In all these systems it is now well established that the environmental cou- 
pling to tunneling system, makes a tunneling motion less frequent and con- 
sequently reduces the effective tunneling matrix element, sometimes termed 
as Debye- Waller factor. This reduction stems from orthogonal properties of 
many particle wavefunctions with different local potentials and is known as 
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Anderson orthogonality theorem [13] . In our study we examine the ground 
state properties of two state system linearly coupled to a phonon bath which 
provides Ohmic dissipation. This system can be described by a spin-boson 
Hamiltonian. 

H = -A a x + u k a\a k + ^ 9k(a k + a ] _ k )a z , (1) 

k k 

where A Q represents bare tunneling matrix element, cr x ,a y are Pauli matri- 
ces. The bath is described by a set of harmonic oscillators and with coupling 
constant g k and a k , a k are the annihilation and creation operators repre- 
senting phonon modes. We have taken coupling constant g k of the form 
gk — 9o | k I 1 / 2 , which represents coupling of a tunneling system to a one- 
dimensional phonons. The frequency dispersion is taken to be uo k =\ k | 
(acoustic phonons) , where k is the quasimomentum wavevector of phonons. 
The Hamiltonian (|IJ) can also be written in terms of creation (cj )and anni- 
hilation (cj) operators of the tunneling particle, namely 

H = -A D (c t 1 c 2 + c{ci) +^2u k a[a k + ^0fc( c i c i - c\c 2 )(a k + a ] _ k ) , (2) 

k k 

where 1 and 2 represent two equivalent positions of tunneling particle. Since 
we are dealing with a single particle it does not matter whether it is a boson or 
a fermion. The Hamiltonian for two state system is also a good representation 
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for statics and dynamics of a quantum particle moving in symmetric double 
well potential ,i.e, in the limit of low temperatures, /cgT much smaller than 
Tiuj {ujq the frequency of small amplitude oscillations in the wells ) and V Q 
the barrier height between the wells [1] . The tunneling frequency A is of the 
order of [\Ju V /h]exp(— 2V /huj ) , in the WKB limit. The Ohmic dissipative 
bath corresponds to a special case [1] where the spectral density S(a;)is given 
by 

S(u)='%2gZ6(u-u k ) = ]-au , (3) 
k z 

where a is the dimensionless dissipation parameter. 

The above Hamiltonian exhibits a rich phenomena for a tunneling system 
ranging from damped oscillations(coherence) up to localization [1] . It has 
been shown that for the weak coupling the particle tunnels back and forth 
between two states and quickly looses the memory of the initial quantum 
state. As the coupling exceeds the critical strength , spontaneous symmetry 
breaking occurs , leading to localization where the effective tunneling matrix 
is identically zero [1,14,15] . This symmetry breaking transition occurs only 
in the limit of zero temperature. At finite temperature and for larger coupling 
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strength the tunneling particle hops incoherently from one state to another. 
The zero temperature localization phenomenon is associated with the infrared 
divergence induced by low energy phonons ,i.e, the phonon field keeps the 
particle localized at the expense of generating an infinite number of low- 
energy phonons. The low freqency phonons play a dominant role in the 
statics as well as in the dynamics of the two-state system. The conventional 
weak coupling picture has a rather limited range of validity and the adiabatic 
approximation over estimates the occurence of infrared divergence. 

Using the renormalization-group procedure the derealization to local- 
ization transition has been studied [1]. When the bare tunneling matrix 
element A Q is small compared to the upper phonon cutoff frequecy cu c , the 
effective tunneling matrix element is given by A e ff = A (=^ 2 -) . At 
a — 1 there is a transition from delocalized to localized (A e // = 0) state. 
The renormalization-group results have been obtained mainly based on the 
assumption of a dilute instanton gas or dilute non-interacting blip approx- 
imation [1] . In the later studies on dissipative quantum tunneling , it has 
been shown that this approximation breaks down due to the condensation of 
instanton gas [16] . It has also been pointed out that the non-interacting blip 
approximation leads to a complete neglect of reaction field which becomes 
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important for the long time dynamics in the larger coupling domain [17]. 
Silbey and Harris [18], and Tanaka and Sakurai [19] have independently suc- 
ceeded in rederiving the renormalization-group results by a simple variational 
procedure which amounts to a kind of mean field theory They have taken 
the displacement of phonons in a ground state via variational parameters. It 
should be noted that when a particle tunnels from one state to another it 
has two distinct effects on the phonon wave function. One is that of a dis- 
placement and the other is that of a deformation. High-frequency phonons 
{uj > A Q ) follow the tunneling particle adiabatically. Hence the displaced 
phonons gives the correct result in this domain. However , for low-frequency 
modes and in the strong coupling regime such an approximation is not suf- 
ficient. The low frequency phonons in general are not always in equilibrium 
with the particle motion , and hence the tunneling particle not only moves 
in an adiabatic potential , but also gets influenced by a retarted potential 
arising due to the phonon bath [4,21]. The motion of the particle at time t' 
disturbs the phonon state , which in turn acts on the motion of the particle 
at later time t. For the frequencies u < A G , these non-adiabatic effects 
dominate , which in turn apart from displacement produce strong deforma- 
tion in the ground state phonon wave function. Chen and coworkers have 
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treated these effects via a variational treatment by introducing a squeezing 
parameter for the phonon states [20,22-24]. Within this treatment they have 
analysed the symmetry breaking transition and have also shown that the 
ground state properties of the tunneling particle not only depends on the 
spectral density , but also on the explicit form of the coupling strength on 
frequency of phonon modes [20]. This non-universal behaviour is at vari- 
ance with the earlier results based on the renormalization-group and the 
obtained phase diagram is in no simple way related to the known results [1]. 
In Ref[25] , the phase diagram for the delocalization-localization transition 
has been analysed through a variational treatment which incorporates both 
the displacement and the squeezing of wave function via two independent 
variational parameters for each phonon mode. 

Most of the theoretical treatments to date are concerned with the effect 
of phonon bath on the tunneling motion. Very few treatments exist which 
directly address the question of effect of tunneling particle on the phonon 
bath [20,26,27]. It has been shown that [20] two state coupled system at 
zero temprature can be mapped on to the modified quantum sine-Gordon 
model , which describes the effective non-linear interaction between the low 
frequency phonons when the tunneling system is in its ground state. It is also 
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of interest to know that the ground state properties of (1 + 1) dimensional 
quantum sine-Gordon model can be mapped on to several two dimensional 
statistical mechanics models ,e.g., two dimensional Coulomb gas etc [27]. In 
the critical regime , where the particle is still in the delocalized state , a gap 
in the elementary excitation of the low frequency phonons is opened near zero 
momentum , and the corresponding ground state wave function is found to 
be of a pairing quasi-particle state , analogous to the BCS super conducting 
state [26,27]. 

In our present work we have carried out a systematic analysis of ground 
state properties of the tunneling particle as well as of the phonons via a 
variational treatment which correctly treats both the displacement and de- 
formation of each phonon wave function in its ground state . As stated earlier 
we have restricted our treatment to a two state system , coupled to bosonic 
enviroment, with a Ohmic spectral density. Results do depend on specific 
form of dependence of coupling constant on the frequency of phonon mode, in 
our case g k = g ^fuj k , this coupling is common to electron-phonon coupling 
in solids. The reduction in the effective tunneling matrix element of the par- 
ticle as a function of dimensionless dissipation parameter a is calculated and 
is compared with earlier known results. The modified phase diagram for the 
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symmetry breaking transtion is presented. We have also shown that in the 
critical regime and in the delocalized phase the gap will be opened near zero 
momentum in the spectrum of the low frequency phonons. The gap vanishes 
throughout in the localized domain. An analytical expression for the phonon 
frequency gap as a function of effective tunneling matrix element and dissi- 
pation parameter a is obtained. We have also obtained an information about 
the modified phonon spectrum and using this , the phonon displacement field 
spatial correlation function has been calculated. In the delocalized phase the 
phonon field correlation functions at large distances (in the asymptotic do- 
main ) falls of exponentially. This implies that in the presence of a tunneling 
particle the lattice is stable under the long wave length fluctuations , so long 
as effective tunneling matrix element for a particle is non-zero. 

2. Ground state variational wave function 

In the following we first motivate our ground state variational wave function 
of the system described by the Hamiltonian given in eqn (1) or (2) . Ground 
state wave function for this system can not be obtained analytically. Only 
special cases can be solved exactly. The simple cases being (i)A a = , and 
(ii) the coupling constant — . In the first case where A D = , the 
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system is just an oscillator and there is no tunneling motion between the 
two degenerate states. However ,the oscillators are displaced in one direction 
when the tunneling system is in one of the two levels and are displaced 
in the opposite direction when the tunneling system is in the other level. 
Correspondingly the total system in its ground state can be represented by 
a doubly degenerate localized state 



Here | 0) represents the vacuum state and i being the position of the tunneling 
state. For an Ohmic dissipative bath the overlap integral between the above 
two states tends to zero , if there is no lower cutoff frequency , i.e, the two 
ground states are orthogonal to each other. In the other cases where gu = 
0,i.e., tunneling particle and phonon modes decouple , the eigenstates of the 
tunneling particle are given by the symmetric and antisymmetric combination 
of two level states , the symmetric combination being the ground state. These 
two situations indicate that for a finite value of and A G the system exhibits 
the competition between the localization resulting from the phonon modes 
and the derealization due to the tunneling. The particle tunnels from one 
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state to another , carrying a dressed phonon cloud. In view of these facts , 
we immediately write down a ground state variational wave function as 

if; = ^Y[exp[-^ | u k | 2 ] exp[-u k al\ exp[-^ k (a k a„ k - a\a}_ k )]C\ 
+ \{exp[~\u k \ 2 } exp[u k a\] exp[--f k (a k a^ k - 4al fc )]C]| | 0) , (5) 



where u k and ^ k are both variational parameters for each phonon mode. We 
call our variational wave function (eqn.[|[]) as VWF(3) . When the param- 
eter 7fc takes a finite non zero value we say that the corresponding phonon 
state (or two phonon state) is squeezed [28]. The squeezed state represents 
a non-classical minimum uncertanity state and are eigenstates of a linear 
combination of the annihilation and creation operators. It differs from a 
minimum uncertainty coherent state in that it has uncertainty less than that 
of a coherent state in one of the quadrature phase given by exp{—^ k ) and in 
another exp^k) , which means that one of the quadrature phase is squeezed 
at the expense of other. In real space representation the displaced squeezed 
state is given as (f)(x) = [irS 2 ]^ 1 ^ exp[— (x— xo) 2 /2S 2 }, where xq is the mean 
displacement and S = exp^k). For this state the uncertainty , Ax in spatial 
coordinate equals to exp{^ k ) and Ap in the momentum equals exp(— 7&). If 
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7fc takes a positive value the wave function spreads over a larger distance 
than the ground state wave function of the harmonic oscillator (coherent 
state). Squeezing parameter takes care of the deformation associated with 
the low frequency phonon wave function as discussed earlier. In VWF(3) if 
we substitute 7^ = , we get the variational wave function proposed inde- 
pendently by Tanaka and Sakurai [19] and Silbey and Harris [18] , we call 
this VWF(l). In the VWF(l) , if we substitue Uk = , it becomes an exact 
ground state of the unperturbed tunneling particle. On the other hand if 
we put Uk = —gk/wk, the wave function describes a linear combination of 
two degenerate ground states given in eqn([|). The VWF(l) does not take 
into account the deformation of low frequency phonons. In VWF(3) if we 
take varitional parameter w fc to be constant and equal to (—gk/^k) , we get 
the variational wave function proposed by Chen et.al [20,22-24] . We call 
this as a VWF(2). The VWF(2) takes in to account the deformation of low 
frequency phonon via squeezing parameter 7^ . However , it does not take 
properly in to account the displacement of phonons. All the phonons are 
assumed to be displaced by an amount | Uk |= gk/^k , which is true only in 
the adiabatic limit or for phonon frequency u > A . In the VWF(3) both Uk 
and 7^ being the variational parameter , naturally it takes care of non adia- 
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batic effects in both the displacement and the deformation of low frequency 
phonons properly. 

Now using VWF(3) one can obtain the expectation value of the total 
ground state energy (E g ) . Using the well known properties [20,22-30] of 
squeezed states and bose operators we get, 

E g (u k ,7k) = -A |exp^[(-2Mfc)exp(-47 fc )]| 

+ u k [u 2 k + (sinh2j k ) 2 ] + ^ 2g k u k . (6) 
k k 

The variational parameters u k and 7^ are determined by minimising E g with 

respect to u k and 7^ , i.e. , we set = and = 0. From this we obtain 

u k = [-g k /(2A Zexp(-A'j k )+uj k )} , (7) 

lk = (l/8)Zn [l + 8A Zu 2 k /u k ] . (8) 

The factor multiplying A D in the eqn. (6) is identified as a tunneling reduction 
factor Z (or Debye- Waller factor ), and is given by 

Z = exp|^[-2^exp(-4 7fc )]| . (9) 

The effective matrix element A e ff is given by A e jj = A Z . The reduction 
factor Z can be determined by using eqns.@ and @ self-consistently. One 
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can notice from eqns.(^) and (||) that u k and ^ k tends to {—g k /uj k ) an d 
zero respectively in the high frequency limit as discussed earlier. This is 
consistent with an adiabatic treatment in this frequency domain , i.e. , high 
frequency phonons with u greater than the A e jj follow the tunneling particle 
adiabatically. 

3. Effective phonon Hamiltonian 

We will now derive an effective Hamiltonian for phonon bath at zero 
temperature consistent with the VWF(3). To this end we first perform a 
canonical transformation 



Ui = exp 
on Hamiltonian ([!]). we get 



(10) 



H 1 = U\HUx = ^uj k a\a k - A a x Cosh ^2u k (a k - a ] _ k ) 
k Ik 

+iA a y Sinh ' s ^2u k {a k — a)_ k ) 

. k 

+ J2(dk + u k u k )a z + 2 ^2(u 2 k u k + g k u k ) 



From the Hamiltonian in eqn.(|ll|) we see that there are two distinct types 
of influence on the phonon bath due to coupling with two-state system. The 
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first is the diagonal term contaning a x , giving static influence due to two level 
system being its ground state (a x = +1). The nondiagonal interaction terms 
contaning a y and a z represents the dynamic influence due to the transition 
in two state system between its ground state and excited state on the phonon 
system. These nondiagonal terms are important only at finite temperatures. 
Hence at zero temperature we neglect these terms along with the constant 
factor appearing in eqn. (|TT|) , thus we have 



ear interactions between the phonons. Here nonlinear interactions are not 
only between the phonon in the same mode , but also in the different modes. 
By appropriately defining the field operators one can readily map the H e ff 
on to modified quantum sine-Gordon model (see ref[20]) in 1+1 space-time 
dimension. It is well known that the quantum sine-Gordon model exihibits a 
Coleman phase trasition as a function of the coupling constant. In Ref[26,27] 
critical conditions for this phase trasition have been derived. It is also shown 
that , in the critical regime the gap will be opened in the excitation spectrum 
near zero momentum . The Ref[27], also elaborates on the connection with 




(12) 



Equation. (^) represents an effective phonon Hamiltonian , having a nonlin 
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Kosterlitz-Thouless classical theory of two dimensional Couloumb gas. In 
these references [20,26,27] static properties of two state systems have been 
studied , by mapping to sine-Gordon model. The results obtained are in 
spirit of VWF(2). However , our results differ in that we closely follow the 
spirit of VWF(3) , which itself is a better ground state variational wave 
function. For the effective phonon system (eqn. ||12|| ) , we apply the second 
unitary transformation 

U 2 = exp J2^(a k a- k - a\a ] _ k ) . (13) 

. k 

This transforms the effective Hamiltonian ( |I2] ) into 

Heff = ^2^k[b k b k cosh4j k + (sinh2^ k ) 2 

k 

+ (l/2)(b{bl k + b k b^ k )(sinh4 7k )} 
-A cosh[J22u k (b k -bl k )exp(-2j k )] , (14) 

k 

where b k and b k are 

b k = Ula k U2 = a k cosh2'j k — a , _ k sinh2^ k , (15) 

b\ = Ula\U2 = a k cosh2 r Y k — a^ k sinh2^ k . (16) 

The transformation (|T3| ) some times refered to as a Bogolibov transformation, 
basically it resacles the phonon normal-mode coordinate and momenta. We 



can expand the hyperbolic function in eqn.(|l4|) by normal ordering phonon 
operators. We choose 7^ such that it minimises the ground state energy as 
given by eqn.(|]). With this choice of 7^ we get , after the straight-forward 
algebra (see also ref.[27]) 

Heff = Y.^exp(A lk )b\b k + E - A o Z[(z0+) 4 /24O - (^ + ) 3 (z</>_)/6 + ....] , 

k 

(17) 

where <p + and 0_ are 

<P+ = J2 u k ex P( 2 lk)bl k , (18) 
k 

4>- = ^u k exp(-2j k )b k . (19) 

k 

The ground states of ( |TTD are displaced squeezed states in the original phonon 
bases. 

cj) = exp{J2lk(a k a- k - 4 a -fc)} exp{-a z ^ u k {a ] _ k - a k )} \ 0) . (20) 
k k 

We see that Hamiltonian ( |17D corresponds to a free field model if the 
tunneling reduction factor Z approaches zero. The effective tunneling ma- 
trix element approaching zero implies the symmetry breaking transition for 
a tunneling particle in that the particle goes from the delocalized state to a 
localized state .The effective tunneling matrix element A e // — > also implies 
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the vanishing of nonlinear interaction terms in Hamiltonian ( |T7| ) as well as 
in the Hamiltonian flT2] ) . In eqn. ( |TTD we identify Q(k) = ujkexp(A~ik) as the 
new excitation spectrum of the phonon quasiparticles with their ground state 
wave function given by eqn. (^) . It is clear that in the critical regime where 
Z is small (Z << 1) , the higher order terms in eqn. ([l7|) are small and cor- 
respondingly Q(k) is well defined excitation spectrum of a modified phonon 
bath in the presence of a tunneling particle. In this limit we later show that 
the gap will be opened in the vicinity of the zero momentum in the excita- 
tion spectrum of the phonons and this gap will disappear as Z — > . The 
opening up of the gap in the vicinity of the zero momentum naturally takes 
care of infrared divergence problem arising in a effective tunneling matrix 
element without invoking infrared cutoff prescription. With a modified exci- 
tation spectrum near the critical regime , we will also calculate the phonon 
displacement field correlation functions. 

4. Numerical procedure 

We will first outline the numerical technique used to investigate the be- 
haviour of the reduction factor from the expression for Z (eqn||) . For the 
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case of Ohmic dissipation this can be written as 

InZ = -a duoujexp(-A-i k )/(2A Zexp(-A^ k ) +uo) 2 , (21) 

Jo 

where a = g^/it is a dimensionless dissipation parameter, which characterises 
the coupling between the particle and the surrounding medium and io c is the 
upper phonon frequency cut-off. Throughout we have taken u k =\ k \. To 
evaluate Z we have to use u k and 7^ as given in eqns.(|7|) and © respectively. 
The expression for 7 fc contains u k and vice-versa. Substituting these values 
for u k and 7^ in eqn. (^1|) leads to an infinite (continued fraction) expression. 
In ref [25] , this self-consistent equation has been solved by iterative procedure 
where one terminates the infinite heirarchy to a finite one. This procedure is a 
approximate one. We have developed an exact numerical procedure without 
resorting to any approximations. To this end we define a new quantity Q(k) = 
Uf t exp(4:'y k ) with this definition eqn.(^Tj) reduces to 

InZ = -a f 1 duott{k) / (2& Z + Q(k)) 2 . (22) 

Jo 

Using eqns (7) and (8) one can readily show that Q(k) obeys a forth order 
polynomial equation , namely 
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n 4 + AA Zn 3 + (AA 2 Z 2 - 8A Zg 2 /to k - cu 2 k )n 2 

(Auj 2 k A Z)n - 4:iu 2 k A 2 Z 2 = . (23) 

Here gk — 9o \ ^ k \ 1 ^ 2 - In eqns.(|22|) and (p3|) we have rescaled the energy 
with respect to u c , for example A Q = A /uo c , Q(k) = Q(k)/u c etc. For a 
given value of a{= gl/ir) , and Z this equation gives four solutions. It turns 
out that out of these four solutions two are complex and one is positive and 
another is negative. We take only the physical solution, i.e., positive one as 
it represents the energy of the modified phonon modes. In some parameter 
range where A o /u is large, four solution to Q of eqn.fl2"B|) are all real. It 



turns out that three of them are negative ,hence we again take a positive 
solution. We first take a fixed value of a , and A D and for a value of Z lying 
in an interval 1 < Z < we calculate Q(k) for all u k . Substituting this 
Q(k) for all uj k in the expression (^) we check for self-consistency. We start 
initially with Z=l and go on changing Z in a small step till we reach the 
self-consistent solution. 



22 



5. Results and discussion 

In fig.(l) we have plotted the reduction factor Z as a function of dimen- 
sionless parameter a for a fixed value of A /u c = 0.1. For comparison we 
have also plotted the earlier results obtained from VWF(l) (....) and VWF(2) 
( ). In all these cases the reduction factor decreases as expected as a func- 
tion of a and at a particular value a c there is a symmetry breaking transition 
from delocalized state to localized state, at which Z becomes identically zero. 
For A /u c = 0.1 the VWF(l) gives a c around 1, whereas VWF(2) gives a c 
around 2. In our case we obtain a c = 1.85. We have identified a c such that 
at this value Z becomes less than 10~ 7 . In parameter range of a between 
zero and 0.5 the results obtained from VWF(l) and VWF(3) are almost 
identical. As we approach a c the results deviate. This deviation should be 
identified as a result of inclusion of fluctuations over the mean field results 
obtained earlier [1]. It is also of interest to note that the parameter range 
< a < 0.5 in the limit A /u c « 1 is of most interest for the quantum 
coherence problem [14]. In this regime at temperature T=0, the tunneling 
system exhibits damped oscillations (coherence) with frequency A e //. The 
symmetry breaking transition obtained from VWF(3) is continous in the 
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sense that the reduction factor Z apporaches zero continously as a tends to 
the critical value a c . Even for larger values of A /oj c we have observed that 
the results obtained for the reduction factor Z from VWF(l) and VWF(3) 
are identical in the range of small dissipation parameter a << 1. In the same 
parameter range of a the results obtained from VWF(2) and VWF(3) differ 
considerably. Moreover, it should be noted that the results for Z obtained 
from VWF(l) in the limit a « 1 are cosistent with the known perturbation 
results [17]. 

In fig(2) we have plotted the phase boundary between the delocalized 
(Z^O) and localized (Z = 0) phase for various values of A q /uj c . In the limit 
of A /u c — > the phase boundary is located at a c — 1. The result based 
on instanton calculation or path integral calculation also predict the phase 
boundary to be at a c = 1 as A Q /u c — > 0. It also predicts the phase boundary 
to be at a c — 1 till A /oj c reaches a value of l/2e = 0.184. In our case 
the phase boundary shifts to the right side (increase in a c ) continuously as 
we vary A Q /u c . The VWF(2) predicts the phase boundary starting from 
a c = 2 in the limit A /u c — > to a higher values of a c as A Q /u c increases. 
In a sense the phase boundary of the derealization to localization transition 
obtained from VWF(l) and VWF(2) are not related. We can say with a 

24 



certan confidence that our results clearly shows the improvement over the 
earlier results obtained by various treatment using path integral formulation 
[1]. In effect our results amount to taking care of flucutations (via a squeezing 
parameter) in an otherwise mean field results obtained earlier. 



From eqn.(|17|) we have identified Q(k) = \ k \ exp{A r )k) as excitation 
spectrum of the modified phonon bath under the coupling with a tunneling 
particle. This excitation spectrum aquires a well defined meaning only when 
the A e jj is small , in this regime the higher order correction terms in eqn. (|T7|) 
becomes quite small and give a small correction to the excitation spectrum 
via perturbation. In fig. (3) we have plotted the excitation spectrum fl(k) 
versus wave vector k , for A /u c = 0.1 and a = (.9599) for which Z=.1005 . 
We have calculated the excitation spectrum with the help of eqn.([23"l) for a 
given value of a and corresponding value of Z which has been calculated by 
self consistent procedure as mentioned earlier. The eqn. (|2~3"D for Q(k) gives 
four solutions and out of these there is only one physically acceptable solution 
which gives a positive value of Q(k). We see from the fig. (3) that the gap 
opening in the vicinity of zero momentum. For comparision we have also 
plotted the excitation spectrum u k =\ k | of the unperturbed free phonon 
modes (in the absence of coupling to a tunneling system). It is clear from 
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the figure that the low frequency phonon modes are affected due to coupling 
as compared to high frequency modes as expected. In fig. (4) we have plotted 
the zero momentum excitation gap O(0) versus the dimensionless dissipative 
cinstant a , for a fixed value of A /oj c = 0.1. The gap has initially a value 
zero up to some a , then it increases and after exhibiting a maxima the 
gap vanishes identically at the critical value of a c . In the regime where the 
gap O(0) is finite , it is given by O(0) = -2A Q Z + 2^f(2A Zair) . Even 
though in fig. (4) we have plotted O(0) versus a for the full range of a where 
tunneling system undergoes derealization to localization transition , the 
gap O(0) is meaningful quantity only in the critical regime , i.e , a close 
to a c (Z << 1). If one would have carried out the calculations in the spirit 
of VWF(2), i.e, by setting the parameter u k = (—gk/^k) in the canonical 
transformation given in eqn.(10), we can readily obtain [26,27] an analytical 
expression for the excitation spectrum of the modified phonon bath. This is 
given by Q(k) —\k\ [1 + (4iraA Z /k 2 )] 1 ^ 2 , with a gap in the vicinity of zero 
momentum O(0) = \fA-nctA Z . The expression for Q(k) depends explicitly 
on a, Z, A , and k. As the results derived from VWF(2) and VWF(3) for Z 
and consequently the symmmetry breaking transition as a function of a differ 
qualitatively, we have not made the comparision between 0(/c) obtained in 
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the spirit of VWF(2) and VWF(3). 

Once we have the full information about the spectrum one can easily 
compute the phonon displacement field spatial correlation function C(x). 
This correlation function is defined as [27] 



where k c is the upper phonon frequency cutoff k c = u c In fig. (5) we have 
plotted C(x) versus x (=xk c ) for A q /lu c = 0.1, a = (.9599) and Z = (.1005). 
The correlation function intially follow a logrithmic behaviour between the 
two space point with the increasing of their space distance x.At the large 
distances (or asymptotically) C(x) decays exponentially after showing a os- 
cillatory behaviour in the intermediate regime. The charcteristic decay length 
£ (or the correlation length) is inversely proportional to Q(0) and hence £ 
becomes large in the critical regime and diverges at a c , the transition point 
of derealization to localization transition. We have also observed on the 
expanded scale that the domain of x « £, the correlation function decays 
logarithmically. It is important to note that in the absence of coupling of 
tunneling particle to lattice phonons , the phonon field correlation function 




(24) 
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diverges logarithmically. This indicates the instability of lattice against long 
wave length fluctuations. However ,by coupling lattice to a tunneling system, 
the lattice acquires stability against the long wave length fluctuations so long 
as tunneling system is in a delocalized state ( i.e , when the excitation spec- 
trum has a gap in the vicinity of a zero momentum ). This fact manifests 
itself in exponential decay of C(x) in the asymptotic regime . Thus a single 
impurity ( or a tunneling system ) can stablize a one dimensional lattice. 
The correlation length £ scales as £ ~ Z^ 1 / 2 and it diverges as we approach 
a c from below. 

In conclusion we have studied the ground state properties of a coupled 
spin-boson system i.e., a tunneling system coupled to an Ohmic dissipative 
interaction arising due to a linear interaction of tunneling system to a one- 
dimensional set of harmonic phonons with an appropriate spectral density. 
We have presented a modified phase diagram for symmetry breaking tran- 
sition. The results for the behaviour of reduction factor on the dissipation 
parameter are given and compared with known results. We have also shown 
that the effective phonon system shows a gap in the excitation spectrum in 
the critical regime and phonon displacement field correlation function decays 
exponentially in the same regime , this stabilizes lattice against the long wave 
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length fluctuations. Throughout our analysis we have restricted to a special 
form of coupling constant , namely is proportional to | k \ 1 / 2 . It is known 
that the different forms for coupling constant lead to qualitatively different 
results for symmetry breaking trasition , in a sense there is no universality in 
a dissipative two state system [21]. The zero temperature properties depend 
not only on the phonon spectral density , but also on the explicit dependence 
of coupling strength on the phonon frequency. We are at present exploring 
the behaviour of effective phonon spectrum for different dependence of cou- 
pling constant and the result will be presented elsewhere. 
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FIGURE CAPTIONS 

Fig.(l). Plot of reduction factor versus dimensionless dissipation constant 

a , for a fixed value A /u c = 0.1. The , , and - - lines are obtained 

from VWF(l), VWF(2) and VWF(3) respectively. 

Fig. (2). Phase diagram for the derealization- localization in (A /u c ) — a 
plane. 

Fig. (3). The elementary excitation spectrum (Q(k)/u c ) of the low fre- 
quency phonons versus scaled wave vector k/k c for a fixed value of A G = 0.1, 

a = .9599 and the corresponding value of Z=0.1005. The dotted lines( ) 

represent excitation spectrum ic^ —\ k |, for uncoupled phonon modes. 

Fig. (4). The zero wave vector gap Q(0)/u c in the excitation spectrum 
versus dimensionless dissipation constant a . For a fixed value of A G = 0.1. 

Fig. (5). The phonon field displacement spatial correlation function C(x) 
versus scaled distance k c x for a fixed value of A G = 0.1, a — .9599 and the 
corresponding Z=0.1005. 
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